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Summary 

An experimental study was conducted to analyze the 
response and evolution of a two-dimensional incom- 
pressible turbulent boundary layer after being strongly 
perturbed by a one-sided sudden expansion. The hot-wire 
measurement technique was used to measure three 
Reynolds stresses and higher-order mean products of 
velocity fluctuations. The Reynolds number, Re h , based on 
the step height, h y and the reference velocity, U 0 , was 
37,200. The upstream flow was a fully developed turbulent 
boundary layer with momentum-thickness Reynolds 
number Re e = 3600 and 5 99 //i = 0.82 . 

The results show that an internal shear layer, imbedded in 
the inner part of the original boundary layer, emanates 
immediately from the step edge. This shear layer has many 
similarities with a plane mixing layer, but does not 
resemble it exactly. 

The new no-slip and impermeability boundary conditions 
imposed on the flow in the reattachment region limit fur- 
ther growth of the mixing-layer-like flow, promoting a 
new internal boundary layer developing on the wall down- 
stream of reattachment. The data show that the structure of 
the internal layer attains quasi-equilibrium, with produc- 
tion and dissipation of turbulent kinetic energy approxi- 
mately equal, by about 20 h from the step. The skin friction 
coefficient has nearly reached its maximum value at this 
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location. In the external layer, the Reynolds stresses, normal- 
ized by the local u \ , initially decrease at about the same rate 
as those in the internal boundary layer. However, the rate of 
recovery reduces sharply beyond x ~ 20/z so that even at the 
last measuring station, x = 5 1 h , the external layer still has 
a memory of the upstream disturbance. 

Nomenclature 

Cf o skin friction coefficient, 2 T^pU 0 

C fe local skin friction coefficient, 2 x w /pU~ e 

C p pressure coefficient, 2(p-p 0 )/pU~ o 

h step height 

X r mean reattachment length 

p static pressure 

Re h step height Reynolds number, U 0 h/v 

Reg momentum thickness Reynolds number, 

U 0 e/v 

Uq upstream freestream reference velocity 

U mean velocity in streamwise direction 

u x friction velocity 



M, V 

fluctuating velocity components in x and y 
directions respectively 

w rms’ v rms 

root mean-square of u and v fluctuations 


normal Reynolds stresses 

~uv 

Reynolds shear stress 


U 2 V, uv 2 

triple-product correlations of fluctuating 


velocity components 

x.y 

coordinate system representing stream- 
wise and wall-perpendicular direc- 
tions measured from the step and the 
wall respectively 

y + 

normalized distance from the wall, yu/v 

V 

molecular kinematic viscosity of air, nom- 
inally 1.5*10“ 5 m 2 /s at T = 20°C 

p 

air density, 1 .2 kg/m 3 at T = 20°C 

5 

boundary layer thickness where 
U = 0.99U e 

So 

boundary layer thickness upstream of the 
step 

8 

8* 

displacement thickness, j* ( 1 - U/U € ) dy 

e 

0 

momentum thickness, 

8 

\{U/U e ) ( \~u/u e )dy 
0 


wall shear stress 


Introduction 

Separating/reattaching flows occur in a wide variety of 
practical engineering applications such as airfoils at angle 


of attack, re-entry vehicles, diffusers, turbomachines, 
combustors, sudden area changes in pipes or ducts and 
atmospheric flows over fences and hills. Separation 
generally has a negative impact on the performance of 
these systems, such as reduced maneuverability of aircraft 
and increased drag, increased local wall heat transfer, 
increase of dynamic structural loads and noise, etc. There 
have been many studies on separating/reattaching flows 
during the past four decades. The research has been 
conducted for different geometric configurations, mainly 
two-dimensional. Fundamental features of this class of 
flows have been addressed most frequently for the 
backward-facing step, the normal flat plate with a splitter 
plate, and the blunt flat plate; in all these cases the 
separation line is straight and fixed by the geometry. 

A separating/reattaching flow can be divided into four 
interacting zones (see fig. 1). The zones are: the separated 
free shear layer, the recirculating region under the shear 
layer, the reattachment region and the attached/recovery 
region. Each flow region bears some similarities to well- 
studied flow cases such as mixing layers and boundary 
layers, while the reattachment and recirculating regions are 
unique to separated flows. 

Different researchers have addressed various issues of this 
class of flows. 

The objective of this report is to present a set of detailed 
measurements of the separated flow behind a backward- 
facing step, and to perform a detailed analysis of the data. 
The analysis is confined to the separated free shear layer 
and the reattachment and attached/recovery regions: the 
hot-wire measurement technique does not give useful 
results in or near re versed-flow regions. The controversial 
issue of the similarity of the separated shear layer to a plane 
mixing layer is assessed in some detail with particular 
reference to the presence of a developing internal mixing 
layer. A thorough analysis of this possible similarity has 
not been carried out in other back-step flow experiments: it 
is more likely to occur here than in the highly-disturbed 
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flow downstream of a normal plate. Secondly, the issue of 
the readjustment of such perturbed flow to the wall 
boundary conditions and its recovery to an equilibrium 
turbulent boundary layer is analyzed. This is just one case 
from a series of five experiments on the backward-facing 
step designed to investigate physics of the recovery of a 
separated flow for 6800 < Re h < 37,000 and 0.8 < 5 99 / 
h < 2.0 (Jovic (1995), to be published as a NASA TM). 

Unsteadiness 

Kim, Kline, and Johnston (1980) observed that the flow 
was highly unsteady in the reattachment region of the 
backward-facing step flow, suggesting violent interactions 
of large structures from the upstream separated shear layer 
with the wall. Two basic modes of char-acteristic 
frequencies are found in all the abovementioned flow 
configurations. The higher-frequency mode is associated 
with the usual large scale motions in the shear layer while 
the lower-frequency mode reflects overall separation- 
bubble growth/decay dynamics or shear layer “flapping” as 
it is frequently called in the literature. In the case of a blunt 
plate, Cherry, Hillier, and Latour (1984), using pressure 
and pressure-velocity cross-correlation measurements, 
observed low-frequency unsteadiness in the separation 
bubble. Eaton and Johnston (1982) argued that the 
observed low-frequency motion in the backward-facing 
step flow is likely to be a consequence of an instantaneous 
imbalance between the entrainment rate from the 
recirculation zone and its resupply near the reattachment 
line. The nature and physics of this imbalance were further 
assessed by Pronchick and Kline (1983). They suggested a 
qualitative model of the reattachment process based on 
flow visualization studies. They observed two groups of 
structures in the reattachment region: “overriding” and 
“interacting” eddies, which will be discussed later. In light 
of the flow visualization and the proposed qualitative 
model by Pronchick and Kline (1983), it can be concluded 
that the separation bubble dynamics depends on very 
complex interactions of large eddy structures, formed in 


the upstream free shear layer, with the wall in the 
reattachment region. 

Kiya and Sasaki (1983), using a similar experimental 
technique to that of Cherry et al. (1984), found that the low- 
frequency unsteadiness has fX r /U ' less than about 0.2 
while the wide-band mode with frequency of 0.6 to 0.8 
corresponds to the frequencies of vortical structures seen in 
plane mixing layers. Castro and Haque (1987) detected the 
same two modes of unsteadiness in the separated flow 
behind a normal flat plate with a long central splitter plate. 
However, Ruderich, and Femholz (1986) observed no 
dominant frequencies in their power spectra for the same 
flow configuration, which led them to believe that there 
was no flapping of the reattaching shear layer. In the case 
of a backward-facing step, Eaton and Johnston (1981); and 
Driver, Seegmiller, and Marvin (1983) also detected the 
two characteristic frequencies with approximately the same 
values obtained by Kiya and Sasaki (1983). Importantly, 
Driver et al. (1983) found that the contribution of low- 
frequency unsteadiness to the total turbulent kinetic energy 
is negligibly small. This suggests that turbulence quantities 
in a backward-facing step flow are not likely to be 
contaminated by low-frequency unsteadiness. 

Similarity to a Plane Mixing Layer 

Cherry et al. (1984) found that the (nearly linear) rate of 
growth of vorticity thickness, the characteristic frequencies 
and the spanwise correlation scales for the blunt plate flow 
are all similar to values found for the fully-developed plane 
mixing layer. Castro and Haque (1987) strongly argued 
that the separated shear layer downstream of a normal plate 
is not closely similar to that of a plane mixing layer. They 
found that the normal stresses are consistently higher than 
in the plane layer. In addition, they concluded that the 
effect of the stabilizing curvature of the shear layer on the 
flow is weak, and that the re-entrainment of the recirculated 
fluid (back into the shear layer) represents the dominant 
mechanism for maintenance of the high normal stresses. 
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Three-Dimensionality 

Cherry et al. (1984) concluded that reattachment affects 
small and large scale flow structures in different ways. 
Their spanwise surface -pressure cross-correlation 
measurements at several streamwise locations indicate 
progressive reduction of length scales near the wall which 
is interpreted as a “three-dimensionalizing of the smaller 
scales.” Spanwise velocity correlations (near the shear 
layer edge where local turbulence intensity was 
2.5 percent) reveal that a three-dimensional state has been 
established well before reattachment and that it is not 
significantly affected by the reattachment. Troutt, 
Scheelke, and Norman (1984) looked at the spanwise 
structure of the flow behind the backward-facing step and 
arrived at the same conclusion, that the small scales are 
subjected to a reduction in size near the wall while the 
large-scale spanwise structure persists far downstream of 
reatta-chment, only to be gradually affected by the new 
boundary condition imposed by the bounding wall. In a 
wall bounded mixing layer experiment. Wood and 
Bradshaw (1982) observed changes in turbulence structure 
before the mixing layer “reattached.” They observed 
significant increases of all spanwise scales as the flow 
became affected by the wall, probably because the 
irrotational “backflow” around erupting eddies in the 
intermittent region, required by the continuity equation, 
tends to be constrained to z-wise motion by the v = 0 
condition at the solid surface. This effect can take place 
without significant changes in the (rotational) turbulence 
itself. 

Near- Wall Region within the Bubble 

A number of researchers such as Westphal, Johnston, and 
Eaton (1984); Adams, Johnston, and Eaton (1984); Driver 
and Seegmiller (1985); Castro and Haque (1987); Ruderich 
and Femholz (1986); and Devenport and Sutton (1991) 
have measured skin friction in the reversed-flow regions of 
different separated-flow configurations. Studies by Adams 
et al. (1984); and Devenport and Sutton (1991) showed that 
the near-wall behavior of the boundary layer within the 


recirculating zone is very different from that of a standard 
turbulent boundary layer. Devenport and Sutton showed 
that mean velocity profiles follow a modified version of 
Simpson’s (1983) model. They also argued that streamwise 
turbulence intensity profiles appear to be scaleable on the 
rms friction-velocity fluctuations and not on the mean . 
Using a simple analytic model, they showed that the extra 
streamwise velocity fluctuations are driven by pressure- 
gradient fluctuations imposed by the shear layer above. 
Jovic and Driver (1994), using laser-oil-film 
interferometry, investigated the effect of Reynolds number, 
on the skin-friction coefficient in the recirculating 

region. They found that in the attached part of the flow 

2 

Cjr = 2t vv /p U Q generally decreased slowly with 

increasing Re, as usual, while the maximum negative C j in 
the recirculating region varies as rapidly as 
suggesting predominantly viscous behavior of the thin 
layer under the recirculating region. The same behavior 
was observed by Adams et al. (1984); and Devenport and 
Sutton (1991). 

The abundance of experimental data obtained for the 
abovementioned flow configurations was restricted mostly 
to the separation bubble and the early part of the 
reattachment region. The results lead to a consensus on two 
aspects of the physics of the separated/recirculating region. 
One is the existence of the two characteristic frequencies in 
the separation bubble and the second is the structure of the 
thin layer near the wall in the recirculating region. It 
appears that none of the experiments studied the very 
important physics of the flow recovery in detail. Cutler and 
Johnston (1989) conducted an extensive study of the 
separated turbulent boundary layer recovery downstream 
of a fence, in an adverse pressure gradient tailored to 
achieve boundary layer equilibrium. The length of the 
recovery region was 83 h ( h is the fence height). They 
indicated that an equilibrium boundary layer structure is 
reached only after several hundred fence heights, which is 
even longer than the value of 100/z suggested by Bradshaw 
and Wong (1972) for a backs tep flow in zero pressure 
gradient. 
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Apparatus, Techniques, and Conditions 

The measurements were performed in a tunnel with a 
symmetric three-dimensional 9:1 contraction, a 169cm 
long flow development section with dimensions 20cm x 
42cm, a backward-facing step of the height, h, of 3.8cm and 
a 205cm long recovery section. The flow was tripped at the 
inlet of the development section using 1.6mm diameter 
wire followed by a 1 10mm width of 40 grit emery paper. In 
order to compensate for blockage effects of the flow due to 
the side wall boundary layers, the side walls were diverged 
according to the estimated displacement thickness of the 
side-wall boundary layers. All the measurements were 


made at a reference flow speed, U 0 , of 14.7 m/s measured at 
a station 40mm upstream of the step. The free-stream 
turbulence intensity as determined by the hot-wire 
measurements in the free stream was 0.4 percent. The 
boundary layer was fully turbulent at a reference station 
1.05/i upstream of the step, having a Reynolds number 
based on a momentum thickness, Re e , of 3600 and a shape 
factor, H, of 1.4. The boundary layer thickness, 8 0 - 5 99 , 
at the reference station was 31mm so that 8 g /h - 0.8. 
The aspect ratio (tunnel width/step height) of 1 1 is just 
above the value of 10 recommended by de Brederode and 
Bradshaw (1972) as the minimum to assure two- 
dimensionality of the flow (in the mean) in the central 
region of a tunnel. The expansion ratio was 1.19 and the 
Reynolds number based on step height was 37,000. 

Surface static pressures were measured on the upper and 
the lower (step-side) walls using a standard pressure 
transducer. The skin friction coefficient distribution 
downstream of the step was measured using a laser-oil 
interferometer. This technique allowed unambiguous direct 
measurements of the shear stress, both in the recirculating 
and the reattached regions of the flow. A more detailed 
description of the method and the results obtained is 
presented in Jovic and Driver (1994). 

Mean velocity and turbulence measurements were made 
with normal and X-wire probes driven by constant- 
temperature anemometers (Miller, Shah, and Antonia 
(1987)) made in-house. The sensor filaments were made of 
10 percent Rhodium-Platinum wire 2.5 | im in diameter and 
0.5mm (or 18 wall units in the upstream boundary layer) in 
length for the X-wire probe, and 1 .25 pm in diameter and 
0.3mm (or 11 wall units) in length for the normal-wire 
probe. The spacing between the crossed wires was 0.4mm 
or 15 wall units. The aspect ratio, l/d, of the sensor 
filaments was 240 for both probes. To improve accuracy of 
the measurements in the regions with higher levels of local 
turbulence intensity, the included angle of the crossed 
wires, customarily 90°, was increased to 1 10°. The 
anemometers were operated at overheat ratios of 1 .3 with a 
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frequency response of 25kHz as determined by the square- 
wave test. The normal- wire signal was low-pass filtered at 
10kHz and digitized at 20,000 samples l sec for 30 sec . The 
X-wire signals were low-pass filtered at 6kHz and sampled 
at 12,000 samples/sec for 30 sec . Analog signals were 
digitized using a Tustin A/D converter with 14 bit (plus 
sign) resolution. The probes were calibrated using a static 
calibration procedure and calibration data of each hot-wire 
channel were fitted with a fourth order polynomial. The 
calibration was checked before and after each profile 
measurement. If the hot-wire drift was more than 
±1 percent of the free stream velocity the profile was 
repeated. In this manner, an error due to mean ambient 
temperature variations was minimized. 

Accuracy of the Hot-Wire Measurements 

It is very important to have a good estimate of the accuracy 
of the data obtained with hot-wires, particularly in flows 
with high turbulence intensity (exceeding 20 percent, say) 
when the turbulence measurements obtained with the 
standard hot-wire technique using “cosine law” begin to 
deteriorate. Accuracy of hot-wire measurements are 
generally affected by uncertainties in all components of a 
chain of instruments used in an experiment: pressure- 
transducer calibration, ambient temperature, hot-wire 
calibration, hot-wire probe alignment, sensor angles, 
sensor length (I/d'), sensor separation, heat loss to 
supports, hot-wire drift, and other second-order 
uncertainties. Using a method of Moffat (1988), and 
Yavuzkurt (1984), the uncertainties in the normal stresses 
m 2 and v 2 , and shear stress ~«v due to the first five listed 
causes was calculated using response equations. The 
maximum uncertainties for each Reynolds stress 
component was found to be not larger than ±5 percent. 

The sensor length and the sensor separation are primarily 
responsible for the X-wire probe spatial resolution. In 
regions with high turbulence levels and small turbulence 
scales (near wall regions) X-wire may incur large errors 
(Nagano and Tsuji (1994)). The sensor length, / = 0.5mm 


of the present X-wires probe is sufficiently long in terms of 
the wire diameter (l/d = 200) to minimize the end 
conduction effects resulting in an approximately uniform 
temperature distribution along the wire (Champagne, 
Sleicher, and Wehrmann (1967)). On the other hand, the 
sensor length should be sufficiently small, of the order of 
Kolmogorov length scale, to avoid spatial averaging of 
small scales along the wire. The length l<5L k (L k is the 
Kolmogorov length scales which is about 0.1mm in the 
near wall region) minimizes undesirable spatial averaging 
(Ligrani and Bradshaw (1987); and Browne, Antonia, and 
Shah (1988)). Ligrani and Bradshaw found that the 
maximum error in u 2 due to l/d and diameter, d, in the near 
wall region could be as high as 7 percent. Sensor separation 
A z ^ 4Z^ may introduce an error of 5 percent according to 
Browne et al., Nagano and Tsuji showed that the most 
sensitive component to the sensor separation is v 2 and that 
error in Reynolds stress components is the function of the 
turbulence intensity. In the separated shear layer region 
where local turbulence intensity exceeds 30 percent, the 
total uncertainty in u 2 is estimated to be ±10 percent, 
±15 percent in v 2 , and ±18 percent in -wv, while the 
uncertainties reduce to less than ±8 percent respectively 
in the recovery region where turbulence intensity is still 
high but gradually decreases. 

High levels of turbulence, exceeding nominally 1 0 percent 
to 20 percent, introduce nonlinear effects into the response 
equation of a hot-wire which cannot be neglected as 
opposed to the standard hot-wire technique (Hinze (1975); 
Bruun (1972), and Muller (1982)). Following the method 
described by Muller, an improved data reduction method 
taking into account triple-velocity products was 
introduced. The truncation error of the series expansion of 
the response equation, which is built into the resulting 
Reynolds stresses, is reduced from third to fourth order. 
The correction due to included triple-velocity products is 
not uniform across the shear layer. The corrected and 
uncorrected Reynolds stresses are shown in figure 1(c). 
Maximum corrections of the Reynolds stresses near the 
wall amount to 18 percent, 35 percent, and 15 percent for 
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u 2 , v 2 , and -«v respectively, while in the outer layer the 
maximum corrections reach the levels of 10 percent, 
35 percent, and 20 percent respectively. It appears that the 
v 2 is most sensitive to the nonlinear effects of the X-wire 
response equation. The correction gradually diminishes 
downstream of the reattachment as the flow recovers from 
the separation. Only corrected Reynolds stresses are 
presented in the report. 

Errors due to rectification of the anemometer signal cannot 
be accounted for. This problem occurs roughly for y < h 
and x<\.2X r where the instantaneous velocity vector 
occasionally reverses its direction or falls outside of the 
angle formed by the X-wire sensors. Tutu and Shevray 
(1975) estimated that shear stress incurs error of 28 percent 
for turbulence intensities greater than 30 percent which 
roughly agrees with the presently applied correction. It is 
not attempted to correct triple- velocity products since the 
correction implies knowledge of all fourth-order moments. 

Due to the aforementioned accuracy problems encountered 
in high intensity turbulent flows, the results of the present 
experiment should be used with caution in the separated 
shear layer region (0 < x < 7/i), where local turbulent levels 
exceed 30 percent. Figure 1(b) shows contours of 
J^/U = 0.3, roughly the boundary of quantitative 
accuracy and qualitative usefulness of the hot-wire 
measurements, respectively. It appears that maxima of all 
Reynolds stresses falls in the high-uncertain region (see 
fig. 3(bHd)) This uncertainty of the data, however, does 
not significantly alter the general conclusions about the 
separated shear layer downstream of the step. 

The boundary layers on the top and bottom walls of the 
tunnel merge for x > 50 h, hence the profiles of different 
turbulent quantities are contaminated by interaction of the 
two layers. 

Results 

In the present study, a low-viscosity oil was used to 
visualize the flow pattern in the separated region and to 


determine the mean reattachment length on the bottom wall 
of the wind tunnel. Flow reattachment occurs at about 
x/h = 6.8 on the centerline. It was observed that the 
reattachment line is not a straight line in the spanwise 
direction but curves upstream near the side walls because 
of interaction with the side wall boundary layers. The 
reattachment line was nominally straight over 65 percent of 
the wind-tunnel width. 

Wall Pressure Coefficient and Wall Shear Stress 

Distributions of the wall-pressure coefficient, 
C p = 2 (p - p o ) /pU 2 Q , along the top and bottom walls 
are shown in figure 2(a). Most of the pressure recovery on 
the step-side of the tunnel occurs within 10 h of the step 
while on the top wall it takes about 20 h. The separated 
shear layer is influenced by the strong adverse pressure 
gradient, by the streamline curvature and by the presence of 
a highly turbulent recirculating flow beneath it. Castro and 
Haque (1987) argued that the re-entrainment of the 
recirculated fluid into the shear layer dominates the 
stabilizing curvature influence on the flow. They were 
studying the flow behind a normal plate, but the comment 
should also apply to the backstep flow where the curvature 
is less. After the flow reattaches, the recovering boundary 
layer evolves under zero pressure gradient. The maximum 
pressure coefficient is about 0.18. 

The distribution of skin-friction coefficient, 

Cjr 0 = 2x w /p U 2 , plotted against x/X r is shown in 

figure 2(b). The wall shear stress, x^, was measured 

directly using laser-oil interferometry (LOI) throughout the 

separating/reattaching region. Downstream of the 

reattachment point, the skin friction coefficient, 
2 

Cj = 2x w /pU e (note that the local free-stream velocity 
U e was used for normalization), was also estimated from 
the Clauser chart, by fitting mean velocity profiles to the 
logarithmic law of the wall. Note that Jo vie and Driver 
(1994, 1995) showed that the log-law is violated in the near 
field of reattaching flows. The Ludwieg-Tillmann 
correlation was also used to estimate skin friction. 
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Agreement in Cy between the three different methods is 
good to within 5 percent sufficiently far downstream of 
reattachment, x > 20 h approximately. The discrepancy 
between the Cy distributions estimated by the two latter 
methods and the LOI technique, shown in figure 2(c), 
clearly demonstrate that the Clauser chart and the 
Ludwieg-Tillmann correlation are not appropriate 
techniques to determine Cy in reattached/recovery flows. 

Separated Shear Layer and Reattachment Region 

Mean flow- The sudden change of boundary condition as 
the no-slip and impermeability conditions are abruptly 
removed at the step leads to a sudden acceleration of the 
flow near y = 0, producing an inflection point in the 
mean velocity profile (see fig. 3(a)). The presence of an 
inflection point leads to a Kelvin-Helmholtz instability 
and the actual rollup of spanwise vortices immediately 
downstream of the step. This behavior is clearly docu- 
mented in a flow visualization movie by Pronchick and 
Kline (1983). The rolled-up vortices do not occur across 
the entire separated shear layer but are confined to a thin 
internal layer imbedded in the original boundary layer 
(fig. 1). The streamlines shown in figure 3(e) were 
obtained by integrating inwards from a reference stream- 
line near the boundary-layer edge, the inclination of the 
streamline to the (known) line y = 8 being assumed 
equal to that in the upstream boundary layer (i.e., no 
change in entrainment rate). The reason for this indirect 
approach is that hot-wire measurements in the recircula- 
tion region are not reliable enough to permit integration 
out from y = 0. It is seen that the average radius of curva- 
ture of the streamline \| / = 0.5 (which starts near mid- 
layer in the upstream boundary layer and is not shown in 
fig- 3(e)) is about 60 h over the interval 0 < x/h < 5, after 
which the curvature reverses but has generally smaller val- 
ues. A typical value of 5 /R is therefore somewhat less 
than the 0.03 which Plesniak, Mehta, and Johnston (1994) 
found to produce significant alteration of turbulence in a 


mixing layer, suggesting that curvature effects in the 
present separated flow are not large, except possibly near 
the surface for a short distance near the reattachment 
point. This is, however, much smaller curvature than that 
of the normal plate studied by Castro and Haque. 

Reynolds stresses- The profiles of Reynolds stresses in 
the separated shear layer are shown in figure 3(b)-(d). 
These quantities, and higher-order products of velocity 
components, are shown in a lab-fixed Cartesian ( x , y) 
coordinate system, not streamline coordinates. Effects of 
the introduced perturbation (sudden expansion and change 
of boundary conditions) on the shear layer are obvious. 
Figure 3 shows that all three measured Reynolds stresses 
increase significantly in mid-layer (y = /z), displaying a 
slope discontinuity in their profiles. Above the discontinu- 
ity, for larger y, the flow remains virtually unaltered by the 
increased turbulence production in the internal layer. This 
can be seen from figure 3(f), where the shear stress profiles 
of figure 3(d) are replotted against stream function and the 
gradual outward propagation of the internal mixing layer 
is clearly seen: the shear stress on a given streamline out- 
side the internal layer continues to change at about the 
same slow rate as in the upstream boundary layer. A plot 
of the value of stream function at the outer edge of the 
internal mixing layer shows that the growth rate for x/h > 2 
is approximately double that for x/h < 1 : evidently the 
internal layer can propagate more rapidly once its stress- 
producing eddies grow to a size comparable with those in 
the outer part of the boundary layer. 

Reliable values of stream function cannot be obtained in 
the lower part of the internal mixing layer, but further 
insight into the multi-layer structure of the separated shear 
layer can be obtained from a plot of the u-component 
skewness S u = w 3 /(u 2 ) 3/2 (fig. 3(g)). S u takes large 
values near the free-stream edges of any shear layer 
(negative on the high-velocity side and positive on the low- 
velocity side) and goes through zero in the 
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maximum-intensity region of a mixing layer. The line A-A 
connecting the left-hand zeroes in figure 3(g) therefore 
marks the high-intensity region of the mixing layer, and the 
lines B-B and C-C joining the two sets of minima mark the 
outer edge of the mixing layer, and the edge of the original 
boundary layer, respectively, the definition of “edge” being 
somewhat qualitative. Note that the high-intensity region 
(or the peak in shear stress) moves inwards with respect to 
y (fig. 3(d),(g)) but outwards with respect to \| f (fig. 3(f)). 

The sharp demarcation between the internal and external 
layers appears to indicate that the two layers contain large 
eddy structures with different dynamics, which 
communicate only through the - presumably small-scale - 
mixing at the interface. The turbulent stresses shown in 
figure 3(b}-(d) increase in the downstream direction, 
attaining almost symmetric distributions about the local 
peak of each quantity (the profiles in fig. 3(f) are far from 
symmetrical because 3 \|//3y = U changes rapidly with 
y). If the step wall were removed, the evolving shear layer 
would be expected to attain the self-similar structure of a 
plane mixing layer. In the reattachment region, the 
presence of the wall is felt by the flow one or two step 
heights upstream of the mean reattachment point, roughly 
where all the turbulent stresses reach maxima. The same 
behavior was observed by Wood and Bradshaw (1982) in 
the case of a mixing layer constrained by a solid wall, and 
by Chandrsuda and Bradshaw (1981); Eaton and Johnston 
(1980); and Troutt et al. (1984) among others in the case of 
backward-facing step. 

There is an unresolved question about the reason for this 
rapid destruction of the turbulent energy downstream of 
reattachment acquired in the separated shear layer: note 
that streamline curvature becomes destabilizing in the 
reattachment region. According to Troutt et al. (1984), the 
decay of the Reynolds stresses in the reattachment region 
coincides with the inhibition of vortex pairing due to the 
close proximity of the bottom wall. Pronchick and Kline 
(1983), based on their flow visualization, observed a large 
number of different instantaneous events and divided them 


into the two major categories: (i) “overriding” eddies that 
pass over the reattachment zone mostly unaltered and 
(ii) “interacting” eddies which are significantly altered 
after interaction with the wall. They divided the latter 
group further into the three subgroups: (1) recirculated 
backflow - an eddy is recirculated after suffering major 
distortion, (2) downstream interaction - an eddy is tom in 
two so that one portion convects downstream while the 
other one provides backflow, (3) lifted backflow - part of a 
recirculating flow (eddy) is lifted by another overriding 
eddy or an interacting eddy. All three processes lead to the 
reduction of turbulent length scales. The presence of such 
eddies with different origins in different parts of the 
reattachment region give rise to a “discontinuity of history” 
which results in a reduction of the correlation between 
velocity components - i.e., a reduction in shear stress. 

Comparison with a Plane Mixing Layer 

The separated shear layer is influenced by the strong 
adverse pressure gradient, the short development length, 
the presence of a highly turbulent recirculating flow 
beneath it and a sheared turbulent boundary layer above it, 
and possibly by the streamline curvature. Thus, the 
separated shear layer cannot be expected to resemble a 
plane mixing layer exactly. As indicated in the two sections 
above, it appears that the separated boundary layer initially 
responds to the perturbation only in a thin layer close to the 
source of the perturbation (in this case, the step lip) while 
the rest of the external layer remains unaffected. The 
internal layer, which develops imbedded in the original 
boundary layer, must bear some phenomenological 
similarity to a plane mixing layer. Therefore, similarity of 
the evolving internal layer to a plane mixing layer is 
examined in this section, and we begin by defining suitable 
scales for the comparison. The vorticity thickness of a 
mixing layer is commonly defined as 
A = A U/(dU/dy) max , where A U = U £ - U min and 
JJ € is the shear layer edge velocity. The minimum velocity, 
U min' on low-speed side of the present shear layer is 

not evaluated directly due to the inherent deficiency of the 


9 



hot-wire technique in the reversed flow regions, t/ . is 
obtained indirectly by fitting the measured velocity profile 
to the well established velocity profile in the regular 
mixing layer, given by 0.5 [ 1 + tank (t|) ] , where 
T| = (y - y ) / A with y c representing the location of the 
velocity gradient maximum. This trial-and-error procedure 
was complete when the best fit to a given analytic velocity 
profile was established. 

The parameters A U and A were used as the normalization 
parameters for further assessment of the similarity of the 
two flows. The growth of the shear layer vorticity thickness 
is shown in figure 4(a) where the solid line represents the 
vorticity thickness growth taken from Castro and Haque 
(1987). The quantity A/28 0 (fig. 4(b)) expresses the half- 
thickness of the internal mixing layer as a fraction of the 
original boundary layer thickness. Figure 4(b) shows that 
the internal layer indeed grows within the original 
boundary layer and that it apparently spreads across the 
entire shear layer before the flow reattaches at about 
x/h = 7 . The present data do not show a strictly linear 
increase in A as is the case for a plane mixing layer. The 
initial growth rate of the internal mixing layer appears to be 
higher than that of a plane mixing layer, while in the 
reattachment region the vorticity thickness actually 
decreases. Similar behavior of the vorticity thickness was 
observed by Castro and Haque (1987). Figure 4(c) shows 
(A U)/U . The maximum value of about 1.1 indicates that 
U min = -0. 1 U which is in good agreement with the laser 
measurements of Driver and Seegmiller (1985). It 
decreases quite rapidly near reattachment, evidently 
because of the distortion of the velocity profile by the 
induced pressure gradient. Once this starts to happen, 
further resemblance between the step flow and a mixing 
layer cannot be expected. The mean velocity profiles in 
self-similar coordinates are shown in figure 5(a). It is seen 
that the mean velocity within the internal layer collapses on 
the self-similar velocity profile of the plane mixing layer: 
this has been largely forced by the fitting procedure 
explained above. 


The profiles of the Reynolds stresses are shown in 
figure 5(b)-(d) in similarity coordinates. Initially, the 
normalized stresses decrease from the initial high values 
adjusting from the boundary layer to an internal mixing 
layer structure. Apparently, the Reynolds stress 
distributions fail to recover to those of the self-similar 
profiles of Bell and Mehta (1990). The observed overshoot 
of the stresses (for x > 6 h) can be attributed to the slow 
response of turbulence to the reduction of the mean rate of 
strain across the shear layer. The mean rate-of-strain field 
reduces rapidly in the reattachment region due to the local 
acceleration of the flow in the presence of the wall. The 
profiles of the Reynolds stresses in the separated shear 
layer would eventually coincide with those of a regular 
plane mixing layer providing that the separated flow had a 
sufficient streamwise length for its development. However, 
the u = 0, v = 0 boundary condition at the wall gives birth 
at reattachment to a boundary layer, which starts to interact 
with the internal mixing layer. 


Higher-order fluctuating velocity products such as u 2 v and 
mv 2 (others are omitted for brevity) show a high degree of 
similarity to those in a plane mixing layer (see fig. 6). The 
restriction of the growth of large structures by the presence 
of the bottom wall leads to the reduction of triple velocity 
products close to the wall (a negative lobe of the 
distributions) as observed in the experiment of Wood and 
Bradshaw. Near, and downstream of, reattachment, several 
further processes may affect the triple products: break- 
down of the large structures, flow interaction with 
recirculating parts of the tom structures and re-entrainment 
of the same. The lobes in the outer region reduce at a much 
slower rate suggesting that the large structures in the outer 
layer remain almost unaffected. 


The shear correlation coefficient R uv = , 

which represents a measure of the efficiency of turbulent 
mixing, and v 2 /w 2 , are compared with the plane mixing 
layer data of Bell and Mehta (1990) in figure 7. In the initial 
stages of the separated shear layer development, R 
attains a value of 0.6 in the mid layer, which is significantly 
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larger than in the plane mixing layer. It appears that the 
turbulent mixing near the step is higher in the central 
portion of the shear layer (-1 <T| < 1 ). However, R uv 
approaches 0.5 in the mid layer of the flow for increasing 
x. This shows that the structural parameter of the separated 
shear layer approaches that of the self-similar plane mixing 
layer. On the high-velocity side (q > 0 ) of the separated 
shear layer (in the outer layer), for T| > 1.5 and x < 2h 
(x > 2.55 ), R attains a constant value of 0.45 which is 
a value typical for a zero-pressure-gradient turbulent 
boundary layer. This suggests that the external flow is 
unperturbed for x < 2h and that it retains its similarity to the 
upstream boundary layer. Downstream of this x location, it 
appears that the external flow region changes its boundary- 
layer character and approaches that of a mixing layer due 
to the stronger interaction of the internal and external 
layers. 

A similar argument applies to the anisotropy parameter, 
v 2 /u 2 . The competing mixing-layer like structure appears 
to prevail gradually over that of the outer boundary -layer 
like structure for x > 2h. The data show that v 2 / u 2 profiles 
are somewhat lower than in the self-similar plane mixing 
layer. Close to reattachment, the profiles on the low- 
velocity side resemble that of a boundary layer suggesting 
that fluctuations of the transverse v- velocity component are 
attenuated faster than the u component as the flow 
approaches the wall. On the high-velocity side (for 
r| > 1.5 ) and for x < 2 h, the given ratio is initially about 0.4 
which is a characteristic value for a plane boundary layer. 
Further downstream, the ratio gradually approaches the 
distribution of an ordinary plane mixing layer, indicating 
that the structure of the mixing-layer-like flow prevails 
over that of a boundary layer. 

Eddy viscosity and mixing length, shown in figure 8(a),(b) 
respectively, evolve in a similar fashion to all of the above 
turbulent quantities. They are low at the first few stations, 
but approach the self-similar distribution of a plane mixing 
layer for x > 5h. Reattachment occurs at x = 6.8/z. 


Time-averaged triple products of the velocity fluctuating 
components constitute turbulent diffusion terms in the 
Reynolds stress and turbulent kinetic energy balance equa- 
tions. Transport velocities of turbulent kinetic energy, 
V k = -A' 1 (^m 2 v + v 3 j/K 2 [^u 2 + v 2 J , and shear stress, 

V uv = uv 2 / wv, normalized by A U are shown in figure 9. 
Values of approximation constants, in case when w 
velocity component is not measured, have been usually set 
totfj = K 2 = 0.75. However, the DNS data of the back- 
step flow (Le, Moin, and Kim (1993)) indicate that the 
above values overestimate corresponding terms, and that 
the values K x = 0.6 and K 2 = 0.7 provide a better 
approximation. The data show that the evolving V k and V uv 
nominally resemble those of the plane mixing layer. It 
appears that the fluxes of the turbulent energy and the 
shear stress are slightly larger in the central portion of the 
separated shear layer where V k and V uv exhibit a self- 
similar character. In the reattachment region V uv rapidly 
becomes more negative for T|<-1 as -uv approaches 
zero near the wall.In summary, the data indicate that the 
perturbation introduced by the separation of the original 
turbulent boundary layer (TBL) at the step lip does not 
affect the separated shear layer across the entire thickness 
instantaneously. Instead, an internal mixing layer is 
generated, emanating from the step lip, imbedded in the 
external TBL, spreading in the transverse direction, and 
carrying information about the new boundary condition 
across the flow. The results indicate that the internal 
mixing layer experiences an “adjustment phase” in the 
near field (for x < 2/i), while the external boundary layer 
structure remains unaltered. For larger x, it appears that 
external and internal layers interact more vigorously, 
leading to a dominant mixing-layer-like structure across 
the entire shear layer. The overshoot of Reynolds stresses 
at locations close to the reattachment (in mixing layer 
similarity coordinates), when compared with the plane 
mixing layer profiles, can be attributed to the slow 
response of turbulence to a change of the mean rate of 
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strain. It should be emphasized that it cannot be expected 
that the mixing layer will fully resemble a self-similar 
plane mixing layer due to its restricted developing length, 
presence of the wall, adverse pressure gradient and 
presence of the highly turbulent recirculating region below 
the separated shear layer. However, the above findings 
point to a very strong qualitative similarity of the internal 
layer to a plane mixing layer. 

This thin internal layer has a turbulence structure 
resembling that of a plane mixing layer. Above this is an 
external layer unaffected by the perturbation and strongly 
resembling the original boundary layer. By about x = 2h 
the internal layer has filled the inner layer of the original 
boundary layer and commences strong interaction with the 
outer external layer while the whole separated shear layer 
moves towards the wall. 

Recovery region 

In the recovery region, the mixing-layer-like structure of 
the separated shear layer encounters a solid wall at 
reattachment and begins to recover to a structure 
characteristic of a plane TBL. As in the case of the flow 
downstream of separation, the response of the turbulence 
structure to the imposed new boundary condition is not 
instantaneous across the entire flow but is achieved rather 
gradually in y as well as in x. The results show that an 
internal layer forms downstream of reattachment as a result 
of a sudden imposition of the no- si ip boundary condition. 
Initially, the internal layer is dominated by the external 
layer dynamics carrying the memory of the upstream, 
mixing-layer-like, flow structure. However, the near-wall 
structure within the evolving internal layer recovers to that 
of an equilibrium plane TBL, as shown below, although 
recovery is far from complete at the last test station. Three 
different basic flow structures, namely that of the mixing 
layer, and those of the wall and wake layers of the plane 
TBL, compete in the recovery region downstream of 
reattachment. It is clear that this type of flow deviates 
strongly from an equilibrium turbulent flow structure in the 


“near field,” i.e., the region just downstream of 
reattachment for X r < x < 20 h (X r = 6.8 h). 

Mean Flow 

Profiles of mean {/-component velocity in the recovery 
region, measured with normal and crossed hot wires, are 
shown in wall coordinates in figure 10(a). It is apparent that 
the velocity profiles close to reattachment do not collapse 
on the universal law of the wall. This is consistent with the 
results of Jovic and Driver (1994, 1995). Note that the u % 
used was obtained from the direct measurements of t 

w 

using the laser oil-film interferometry technique. The 
profiles exhibit smaller velocity gradients, suggesting 
larger turbulence length scales as discussed by Bradshaw 
and Wong (1972). The large “wake” contribution to the 
mean velocity profiles is not produced by an adverse 
pressure gradient - pressure gradient is negligible in the 
recovery region as seen in figure 2 - but is generated by the 
mixing-layer-like flow upstream of reattachment. This 
“mixing layer” contribution gradually decays with x as the 
outer flow is increasingly influenced by the spreading of 
the newly developing wall layer. This feature of the flow 
during the recovery phase will be discussed more in the 
next section. Forx/h > 20, the mean velocity profile rises to 
follow the universal log-law over a range of y that 
gradually increases downstream, although even at the last 
measurement position, x = 5\h, the profile still dips below 
the log law; the “wake” component is much smaller than in 
an equilibrium TBL and is still decreasing. 

Reynolds Stresses 

Distributions of the three measured Reynolds stresses in 
wall coordinates are shown in figure 10(b)-(d). All three 
stresses attain maxima in the outer layer of the recovering 
boundary layer downstream of the mean reattachment 
point. The peaks have decreased somewhat from the values 
in the separated shear layer but are still several times larger 
than those in an equilibrium TBL. In the reattachment and 
recovery regions, the flow accelerates near the wall, 
resulting in the gradual reduction of the velocity gradients 
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across the outer part of the flow and in an increase of the 
gradients near the wall. Simultaneously, all the Reynolds 
stresses monotonically decrease across the entire 
recovering boundary layer. 

Measurements very close to the wall revealed that the 
normal stress does not scale on the wall variables 
during the initial stages of the flow recovery, as seen in 
Figure 10(b). Jv 2 (fig. 10(c)) is also nonuniversal, at least 
for yu^/v > 30 , but less strongly, -uv is of course forced 
to collapse just outside the buffer layer, where the total 
shear stress is still closely equal to the wall shear stress, 
nonuniversality of Ju 2 and Jv 2 can be attributed to the 
“splat effect,” i.e., the v = 0 boundary condition 
redistributes fluctuations to u and w components induced 
by the constraint of the large-scale eddies of the former 
mixing layer, for values of y much less than the typical 
wavelengths of those eddies. However this does not 
explain the milder nonuniversality of A The 
w-component turbulence intensity exhibits a plateau, 
u /u =4.5, for 40 < y«_/v < 100 at x = 9.87/z. 
u /u is attenuated for smaller normal distances from 

rms X 

the wall while it increases for larger distances until it 
reaches a maximum in the outer layer. The plateau region 
can be interpreted as a region where the wall and the 
external-layer influences overlap. In other words, the value 
of y at the inflection point is related to the thickness of the 
internal layer. This will be discussed in more detail below, 
with reference to the higher-order statistics, skewness and 
flatness, and the quadrant analysis of the shear stress. It 
appears that the plateau and the associated inflection point 
propagate away from the wall with downstream distance, 
indicating the thickening of the internal layer and the 
outward spreading of the wall influence. 

In the near field of reattachment, the normal stresses decay 
rapidly, while further downstream the rate of decay is 
rather slow. Figure 1 1(a) shows the variation with x of the 
maximum value of u /u in the outer layer and its value 
in the internal layer at yu^/v- 20. It is seen that the 
turbulence intensities in the two layers, the internal 


boundary layer and external layer, are initially 
proportional. At y + = 20, however, ^ rm / u x falls to a value 
of about 3.0 at x = 20/z and remains constant for the rest of 
the recovery region, suggesting that the turbulence 
structure recovers to that of an equilibrium TBL in a very 
thin layer close to the wall, regardless of high turbulence 
levels in the external layer. At the same x location 
(x = 20 h), the turbulence intensity attains a maximum 
value of about 3.5 w x in the outer layer, or in terms of the 
free-stream velocity, U e , and the local skin friction 
coefficient: u rms /U e = 3.5 JcJ/2. We can infer a rough 
general rule that any flow perturbation (in the mid layer) 
equal to or greater than this threshold level will be 
sufficient to alter the near-wall boundary layer structure 
significantly. 

At x = 20 h , all Reynolds stresses attain their 
characteristic equilibrium values (in wall units) near the 
wall. It appears that Cj has nominally recovered 
(fig. 2(b)), and the law of the wall for the mean velocity 
profile has reached as far out as y + = 100 (fig. 10(a)), at the 
same x location. Based on these characteristic features, it 
can be concluded that the internal layer qualifies as a 
quasi-equilibrium layer for x > 20 h. Most likely, the first 
and the second derivatives of the ^ rms ^ u T distribution are 
equal to zero at y + = 20 at this x location. It appears that 
once the quasi-equilibrium internal layer forms, the law of 
the wall holds almost independently of the external flow 
conditions, providing u rms /u x < 3.5 or 

U rms /U e = ^JC/2, where “rms re P resentS 3 
maximum turbulence intensity in the external shear layer. 

The implication is that, in the case of an equilibrium 
turbulent boundary layer with C j = 0.003, a free-stream 
turbulence level of the order of 15 percent, and greater, 
affect the near-wall structure for > ,+ < 20. In the outer 
region, the turbulence intensity continues to decay, as seen 
in figure 11(a), indicating that the turbulence in that region 
has not approached an equilibrium state even at the last 
measuring station, x - 5\h. 


13 



— 2 

The decay of the maximum shear stress, -wv/ U , and the 

2 

maximum turbulent kinetic energy, k/ U e , is plotted on a 

log-log scale in figure 11(b). The peak values of both 

quantities lie approximately along straight lines forx>30h. 

- 1.33 

The peak shear stress decays approximately as x 

though its asymptotic value is nominally unity (not zero) 

while the turbulent energy decays somewhat slower, as 
-1 27 

x ' , suggesting slightly faster shear stress recovery. It 

appears that the rate of energy decay of the present 
experiment is approximately same as that of the 
homogeneous turbulence energy-decay rate for x > 30 h 
(Gibson and Dakos (1993)). 

In summary, two regions are identified, the near-field 
X r < x < 20 h and the far- field x > 20 h . In addition, two 
distinct layers in the transverse direction of the reattached 
boundary layer can be identified downstream of the 
reattachment region. One is very thin, of the order of 20 to 
100v/w T in the near-field, and is identified as an internal 
boundary layer. The internal layer is dominated by the 
external layer dynamics in the near-field region, attaining a 
quasi-equilibrium state by x ~ 20 h . The outer layer is 
characterized by a strong memory of the upstream 
perturbation, i.e., it is characterized by the mixing-layer 
type of Reynolds- stress production and the large “mixing- 
layer contribution” seen in the mean velocity profiles. 
However, for x > 20h the inherited mixing-layer structure 
of the external layer comes under the increasing influence 
of the boundary-layer-like structure of the internal layer. 

Driple Products 

Profiles of u 2 v , v 3 , wv 2 , and w 3 , normalized by f/ 3 , are 
shown at the different measuring stations in figure 12. 
Since we are no longer concerned specifically with the 
internal wall layer, these and later quantities are plotted 
against y/8 instead of Magnitudes of the normalized 
triple velocity products near reattachment are about 
20 times greater than those of an equilibrium TBL. This in 
its own right shows the importance of an accurate model- 
ing of turbulent transport (diffusion) terms in the balance 


equations of Reynolds stresses. Beside high magnitudes, it 
is striking to observe presence of large negative lobes in 
w 2 v and v 3 in the attached region at least as far as 
x = 20 h , despite the attenuating effect of the solid sur- 
face. This shows that the flow still has a very strong mem- 
ory of the upstream mixing-layer-like flow. It appears that 
« 3 is much less sensitive to the newly imposed boundary 
conditions, judging by its slow rate of decay near the wall 
as seen in figure 12(d). The large eddies are affected by the 
proximity of the wall exerted through the u = 0 and v = 0 
boundary condition. The turbulent fluxes associated with 
the diffusion of turbulent kinetic energy normal to the wall 
are severely attenuated close to the wall (y/ 8 < 0.1 say). 

The local minimum and maximum values of the triple 
products occur at approximately 0.28 and 0.68 with the 
zero point at about 0.358: compare the behavior of the 
w-component skewness upstream of reattachment in 
figure 3(g). The reduction of the triple products near the 
wall suggests that the large eddies break down or get flat- 
tened, producing eddies of smaller length scales. This is 
consistent with the flow visualization observations of 
Pronchick and Kline (1983), and with the above findings 
for the normal stress. The negative lobes of w 2 v , v 3 , and 
positive lobes of wv 2 disappear for approximately 
x > 25h , when the triple products resemble those of an 
equilibrium TBL near the wall. The difference between the 
two distances, 20h and 25 h> the former being deduced 
from the normal stress (Ju 2 ) and the latter from the triple 
products, indicate that the contribution of the turbulent 
transport (diffusion term) to the Reynolds stresses for 
x > 20h is not significant. 

I\trbulence Structural Parameters, Higher 
Order Statistics and the Quadrant Analysis 

The structural parameters of the recovering flow are 
compared with those of the upstream undisturbed TBL in 
figures 13-15. Structural parameters v 2 /w 2 and R are 
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shown in figure 13. The low values of both parameters 
near the wall in the initial stages of the recovery region 
indicate relatively stronger attenuation of the v-velocity 
component when compared to the w-velocity component. 
This is consistent with the observations made above for 
the Reynolds stresses and triple velocity products in the 
near field. In the outer layer, values of v 2 /u 2 and R uv 
appear to be generally larger than the distribution 
characteristic for an undisturbed TBL. v 2 /u 2 increases 
while R decreases with increasing x , attaining constant 
values of about 0.5 and 0.42 respectively. The trend of 
increasing v 2 /u 2 suggests somewhat slower recovery of 
the v component when compared to the u component in 
the outer layer. The mildly decreasing trend of R uv 
indicates a faster rate of recovery of the shear stress (or 
reduced efficiency of maintenance of shear stress) 
compared to that of the turbulent kinetic energy (not 
shown here). This is consistent with the rates of decay of 
the shear stress and the turbulent kinetic energy presented 
in figure 1 1 . 

The ratio - uv/v 2 (fig. 13(c)) shows rather clearly how the 
internal boundary layer starting at the reattachment point 
spreads out into the former mixing layer. -uv/v 2 is just 
R uv Ju 2 / v but as a diagnostic of shear-stress behavior it 
has the advantage over R uv that it is not contaminated by 
u and w component fluctuations resulting from the “splat 
effect” of wall constraint mentioned above, which contrib- 
utes to turbulent kinetic energy but not to -uv or v 2 ). The 
value at y /8 = 0.15 changes more slowly than the 
remainder of the profile, because of conflicting effects: 
close to the surface, (the closest data point is at very 
roughly y + = 35 , so the fact that -uv/v 2 ~ 1/y in the vis- 
cous sublayer is immaterial) -uv/v 2 rises slowly from 
nominal zero at reattachment, but even at x/h = 51.18 the 
undisturbed boundary layer trend is approached only for 
y/ 5 = 0.03 approximately. A rising tide of -uv/v 2 
starts to fill up the original profile, whose inner part hap- 


pens to be almost linear on the semi-log plot of 
figure 14(c). When this tide reaches y/ 5 = 0.15 (some- 
where between x/h - 15.13 and x/h = 20.29) the peak in 
-uv/v 2 near y/ 5 = 0.5 starts to decrease. Remarkably, 
the peak is still decreasing strongly at the last station, 
x/h = 51.18 , and is well below the boundary-layer value 
to which it must revert downstream. A possible explana- 
tion of the decrease is that as the internal boundary layer 
grows, it starts to interfere with the large eddies in the 
outer layer, whose structure is qualitatively that of a mix- 
ing layer: evidently, the interaction between two different 
eddy structures produces a less efficient structure than 
either - i.e., it reduces the shear stress for a given turbulent 
intensity. 

The relatively large positive w-component skewness S u 
for y<0.38, shown in figure 14(a), indicate that the 
^-component is lower than its mean value more often than 
not. However, the condition u = 0 implies the presence 
of rare but large positive spikes of ^-component, which 
contribute to the long tail of the velocity probability 
distribution and tend to produce a large flatness factor 
also. There are two types of motions which can be 
associated with positive w-spikes, depending upon the sign 
of the v-component. One is the type of event which occurs 
in quadrant I (Qj) of the (w, v) plane, i.e., u >0 and v > 0, 
which like a Q 3 event is called “interactive” motion, and 
the other is the Q 4 event with u > 0 and v < 0, which is 
usually called a “sweep” (see Willmarth and Lu (1972); 
Wallace, Eckelmann, and Brodkey (1972) for a discussion 
of quadrant analysis). In the early stages of the flow 
recovery, surprisingly, the S v profiles cross zero at two y 
locations, y ~ 0.038 and y ~ 0.38 as seen in figure 14(c). 
For y < 0.038 , positive S y and S u suggest the presence of 
positive u and v spikes associated with quadrant I events 
(interactive motion). This type of motion possibly occurs 
when a sweep (w > 0 , v < 0 ) is deflected by the wall and is 
being ejected back into the flow (w > 0 , v > 0 ). 
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Simultaneous large values of F u and F near the wall 
(see fig. 14(b),(d)) indicate large excursions of the 
fluctuating velocity components from their mean. Note 
that the flatness factor is roughly equal to 3/y , where y is 
the intermittency factor (Corrsin and Kistler (1955)). 
Hence, the large values of F ' and F y in the inner layer 
indicate lower intermittency factor which is produced by 
the large structures which, on the other hand, 
communicate between the outer and inner layer of the 
recovering boundary layer. 

It appears that the S y cross-over point closer to the wall 
coincides with the plateau observed in the streamwise tur- 
bulence intensity u 2 (fig. 10(b)), i.e., it roughly corre- 
sponds to the inflection point near the wall. This first zero- 
crossing of S y gradually moves away from the wall while 
the second crossing (further away from the wall) moves 
towards the wall at the downstream stations, until S y 
becomes positive across the entire boundary layer at about 
x = 28 h. Relatively large negative values of S v for 
0.038 < y < 0.38 for a flow in the early stages of the 
recovery indicate a behavior quite different from a plane 
TBL. Positive S u and negative S v show that Q 4 dominate 
in this y interval. The sweeps intermittently transfer high 
momentum fluid to the inner layer, producing high normal 
stresses, while shear stress remains surprisingly low (see 
fig. 10(d)). Initially, both skewnesses, S u and S , change 
sign at y = 0.38 . This location corresponds to the distance 
where all Reynolds stresses attain their maximum values. 
Ueda and Hinze (1975) found that the location where S u 
change sign coincide with the location where u 2 attains a 
local maximum value in the wall region of an equilibrium 
TBL. Simpson, Chew, and Shivaprasad (1981) found that 
S u crosses zero at y « 0.48 where u 2 has a local peak in 
case of the separating turbulent boundary layer. For 
y > 0.38 , negative S u and positive S y indicate occurrence 
of large negative and large positive u and v peaks, respec- 
tively, suggesting that Q 2 dominates farther away from the 


wall. It appears that the higher order statistics of u- and 
v-components in the inner layer approach those distribu- 
tions of a plane TBL for x > 2 5h coinciding with the dis- 
tance deduced from the turbulent diffusion fluxes. 

Beside the qualitative observations made above for the 
higher order statistics, better understanding of the physics 
of the recovering boundary layer can be gained through a 
quadrant analysis of the measured shear stresses. The frac- 
tional contributions to the Reynolds shear stress from each 
quadrant in the (u, v) plane are determined for two loca- 
tions, x = 9.87 h and 38.55 h, and they are compared with 
the corresponding distributions measured in the upstream 
plane TBL. The results are shown in figure 15. At 
x = 9.87/z (fig. 15(a)) the turbulence near the wall is still in 
a disordered state and the shear stress is still small for 
y < 0.058 (compared to the turbulent energy, say), so that 
contributions from the four quadrants almost cancel and 
the customary normalization by -mv is not very meaning- 
ful. It is more useful to normalize by the Q 2 contribution 
(u < 0, v > 0), i.e., ejections from near the surface, which 
can be regarded as locally-generated turbulence as distinct 
from the Q 4 sweeps imposed from above, and this is done 
in figure 16. Compared with the upstream boundary layer 
there are now stronger sweeps in the inner layer through- 
out the recovery region, and the negative contributions, Q t 
and Q 3 , increase greatly as y decreases. It is curious that 
Q 3 (incoming low-speed fluid) is smaller than Qj (high- 
speed fluid moving away from the surface, which would 
be expected to be a rarer event). Possibly the Q1 contribu- 
tions are the return flow from the “splat effect.” By 
x/h = 38.55 (figs. 15(b) and 16(c)) the quadrant contribu- 
tions look much more like those in the upstream boundary 
layer. 

The fraction of time the uv signal spends in each quadrant 
is also of interest, the relevant quantity being the departure 
from 1/4 (fig. 17). In the upstream boundary layer, sweeps 
from above (Q 4 ) occupy more time than ejections from 
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below (Q 2 ), both near the surface and in the outer intermit- 
tent region, and Ch and Q 3 occupy shorter times except 
near the outer edge. At x/h = 9.87 (fig. 17(b)), Q 4 sweeps, 
though now making a larger contribution to uv than Q 2 in 
the first 30 percent of the shear layer thickness, are shorter 
than the Q 2 events, whereas Q 4 events become much 
longer than Q 2 in the outer part of the flow, occupying 
nearly half the total time. The percentage of time occupied 
by Qj and Q 3 events is roughly the same as in the 
upstream boundary layer. The broad outer-layer peak in 
Q 4 time is still present at x/h- 38.55, where the quadrant 
contributions to uv have nearly returned to the boundary- 
layer distribution. Once more we see the long persistence 
of the mixing layer structure in the outer layer. 


terms. The uncertainty of individual terms is rather diffi- 
cult to estimate; however, the terms containing gradients 
with respect to x are clearly more uncertain on a percent- 
age basis than those containing gradients with respect to y 
though the latter are larger. The uncertainty of the dissipa- 
tion when determined from the difference is rather high. 
Since the spanwise velocity component was not measured, 
the following approximations are introduced. The turbu- 
lent kinetic energy, k , was approximated by 0.7^ u 2 + v 2 J, 
turbulent diffusion in the streamwise direction, uk , and in 
the transverse direction, vk , were approximated by 
0.6^m 3 + mv 2 J and 0.6^m 2 v + v 3 J respectively as sug- 
gested by the DNS data for the separated flow, as indicated 
above. 


Transport of the TYirbulent Kinetic Energy 
and Shear Stress 


The balance of the turbulence kinetic energy at four 
representative streamwise locations (9.87/i, 11.847*, 

20.29 h, and 38.55/i) is shown in figure 18. Distributions at 
six other additional locations are not shown, for brevity, 
but will be used in the discussion below. The turbulent 
kinetic energy equation for two-dimensional flows may be 
written as follows 


= 

ox dy 


-^-u( k + - k + - 

dx V p) dy V p 
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Contribution by fluctuating pressure-velocity covariances 
to the turbulent transport (diffusion term) is typically small 
in wall bounded flows and was therefore neglected. How- 
ever, this approximation may be quite crude in the reat- 
tachment region of separated flows where large pressure 
and velocity fluctuations take place. 


All terms of the transport equation were evaluated from 
the measured turbulent quantities except the rate of dissi- 
pation, which was obtained by difference of all the other 


At reattachment, part of the mixing-layer-like flow is 
deflected upstream (at least in a time-average sense), 
bringing highly-turbulent fluid close to the wall. The tur- 
bulent kinetic energy balance in a true mixing layer has 
very large diffusion terms, the loss by diffusion in the cen- 
tral part of the layer being almost as large as the dissipa- 
tion, with correspondingly large gains near the edges. In 
the reattachment region, the near-equality of gain by diffu- 
sion and loss by dissipation near the wall makes the energy 
balance unlike that in the inner layer of a boundary layer 
(production = dissipation) or in any part of a mixing layer. 
It is seen that the contributions to the turbulent energy by 
the production and convection are negligibly small in this 
part of the flow. Although the large eddies which carry out 
the turbulent diffusion are distorted by the v = 0 boundary 
condition at the surface, gain by diffusion is a major term 
in the energy balance near the wall as far downstream as 
x/h ~ 15 (fig. 18(a) and (b», and there is still significant 
gain by diffusion near the wall at x/h = 38.55. It takes a 
very long time for the flow to regain energy equilibrium 
(production approximately equal to dissipation, as it is 
over most of the thickness of a constant-pressure TBL). 
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It was found that contributions from the three terms, 
longitudinal turbulent diffusion, duk/d jc, production by 
the normal stresses, -^w 2 - v 2 'jdU/dx, and the mean 
flow transport are significant downstream of reattachment, 
in contrast to a regular TBL where contributions by these 
three agents are negligibly small. Production by the 
normal stresses near the wall is negative, since dU/dx is 
positive due to the flow acceleration and u 2 is larger than 
v 2 . Note that advection is small - rapid changes in mean 
flow do not necessarily mean rapid changes in turbulence. 
This negative production of turbulent energy is nearly 
equal to the shear production -uvdU/dy , in the inner 
layer (which is still of the order of magnitude of u^/ (Ky) 
as predicted by the law of the wall) so net production is 
much smaller than dissipation in the inner layer, at least as 
far as x/h = 12 (fig. 18(b)). These features lead to the 
conclusion that turbulence models which use equilibrium 
assumptions are not applicable in the near field of 
reattachment. Note that the ^-equation terms shown in 
figure 18 were made dimensionless by f/^/8 where U e is 
a boundary layer edge velocity and 6 is a local boundary 
layer thickness. 

The distributions of all the terms of the ^-equation in the 
outer portion of the flow strongly resemble the distribu- 
tions of the same terms in an ordinary plane mixing layer 
(fig. 18). Even though the magnitude of each term reduces 
in the downstream direction, this similarity is retained 
even up to about 20 h. This shows that the mixing-layer- 
like structure is still present, but decays gradually in the 
outer part of the flow. On the other hand, the competing 
wake-like structure of a regular TBL asserts its presence 
through the boundary condition u = 0 and v = 0 at the 
wall. It affects the mean-velocity gradients, production of 
turbulence and the turbulent transport through the confine- 
ment and break up of large eddies produced in the sepa- 
rated region. As a result, the peak of the production shifts 
from the central region of the layer to the wall region 


where the velocity gradient begins to dominate and the 
shear stress gradually increases. Transport by the turbulent 
diffusion gradually decreases in the outer region (see 
fig. 18(b),(c)). while the maximum loss moves towards the 
wall. The peak of turbulent diffusion in the central region 
of the flow occurs approximately at 0.56 . The peak of the 
diffusion is surprisingly large, almost as large as the sec- 
ondary peak in production even at x/h = 38.55. Apparently, 
large eddies transfer turbulent energy from the central, 
energy-rich, region toward the wall and outwards to the 
boundary layer edge as can be seen in figure 18(a),(b). The 
data suggest that an accurate prediction of the separated/ 
reattached flow will require a good model for the turbulent 
diffusion terms. 


The production peak occurs at about 0.456 in the trans- 
verse direction and is larger than the near-wall production 
in the near field (up to 20 h). Downstream of this location, 
the production in the wall region rapidly increases. It 
appears that the competing outer layer/mixing-layer like 
structure ceases to affect wall structure by about x = 20 h 
when the familiar wall mechanisms prevail in the produc- 
tion of the wall turbulence. By the stream wise distance of 
about x = 30/z, the transport terms of the turbulent kinetic 
energy further decrease in the outer part of the flow as seen 
in figure 18(b)-{d), gradually approaching the structure of 
an equilibrium constant-pressure TBL. 


The balance of the shear stress, -uv , in four stream wise 
locations (9.87/i, 11.84/*, 20.29/*, and 38.55 h) is shown in 
figure 19, again with the terms made dimensionless by 
f/ 2 /6. The shear stress transport equation for two- 
dimensional flows may be written as follows 





pVdy dxJ 
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The turbulent diffusion due to the pressure- velocity corre- 
lations, d (vp/p) /dx and d {up/p) /dy, and the pro- 
duction by the normal stress, u 2 (dV/dx ), were 
neglected. It appears that the advection term is smaller 
than for the ^-equation. The three terms which dominate 
the transport of shear stress are production, v 2 { dU/dy ) , 
combined longitudinal and transverse turbulent diffusion, 
and the pressure-strain “redistribution” term. In the central 
portion of the recovering flow, the maximum shear-stress 
production occurs at the same location as the maximum 
production of turbulent energy. In the near-wall region, the 
gain by diffusion and the loss by the pressure-strain term 
are not balanced as would be expected from the behavior 
of the ^-equation near the wall. The shear-stress diffusion 
term is of the same magnitude and sign (gain) as the pro- 
duction term in the wall region. The production of the 
shear stress has significantly increased near the wall by 
x = 20 h (see fig. 19(b)). Distributions of diffusion and 
pressure-strain terms in the wall region resemble those in a 
plane TBL by about x = 30 h, which is consistent with the 
behavior of the ^-equation. In the outer part of the flow, all 
three terms decrease significantly beyond 30h so that by 
x= 38.55/i (see fig. 19(c)) the distribution of all terms 
strongly resembles that of a regular TBL. The faster recov- 
ery towards local equilibrium is partly illusory: the “rapid” 
part of the pressure-strain term acts in direct opposition to 
the production term, and is often so modelled: if the two 
terms are merged to give a net production, the shear-stress 
balance looks much more like the energy balance, with the 
“slow” part of the pressure-strain term playing the same 
part as energy dissipation. 

Derived Quantities and their Implications for 
Modeling 

In the recovery region, the shear stress and mean 
streamwise velocity, measured in the Cartesian coordinate 
system, were used to evaluate PrandtPs mixing length, 
l - J-uv/ (dU/dy) , and the eddy viscosity. 


v t = -wv/ (dU/dy) . Though lacking physical 
meaning, these two quantities have been successfully used 
as the basis of data correlations in turbulence models for 
calculating slowly-evolving flows. However, such 
correlations fail in more complex flow configurations, 
such as separating/reattaching flows, where the Reynolds 
stresses respond slowly to the rapid changes of the mean 
rate of strain, nondimensional mixing-length, // 6, and 
eddy- viscosity, v { /(U e S j) (where 8j is the local 
displacement thickness), are shown in figure 20(a), (b) 
respectively, compared with the distributions in the 
upstream fully developed TBL. The values of the 
Reynolds number Re e in the recovery region for the 
given Re h are typically over 9000 compared to 3600 in 
the upstream boundary layer, but mixing length and eddy 
viscosity are not expected to depend significantly on 
Reynolds number in this range. Near the wall the mixing 
length is a linear function of the normal distance from the 
wall. However, the slope is significantly larger than the 
value of k = 0.41 typical for a zero-pressure-gradient 
equilibrium boundary layer. At the initial stations 
( x/h < 20 ) the mixing length reaches a maximum value at 
y/8 » 0.2 (see fig. 20(a)) and then drops sharply, rather 
than becoming nearly independent of y as in a regular 
TBL. The drop suggests the presence of two distinct 
structures with different origins, one being associated with 
the developing internal TBL and the second one being 
associated with the inherited mixing-layer-like structure in 
the external part of the flow. 

The mixing-length slope in the linear inner region gradu- 
ally decreases with downstream distance (fig. 21(a)), while 
the nondimensional mixing length in the outer part of the 
flow rises well above the value of 0.085 characteristic of a 
zero-pressure gradient equilibrium TBL remaining 
approximately constant at about 0.20 in the outer region 
for x > 30 h. Similarly, the eddy viscosity, v f / {U e 5j) , 
deviates both in the wall and the outer flow regions from 
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the distribution of the upstream TBL (see fig. 20(c)). It 
appears that, as in the case of the mixing length, the slope 
near the wall is larger than in the upstream TBL. The slope 
gradually decreases downstream approaching the value of 
kw x . Maximum values of the normalized eddy viscosity 
are shown in figure 21(b). It is seen that 
the nondimensional eddy viscosity in the outer part of the 
flow rises initially, reaching the value of 0.075 at approxi- 
mately x = 20 h. Subsequently, it begins to decay at a very 
slow rate. It is seen that even at the last measuring station, 
x = 5lh y the maximum eddy viscosity is about four times 
larger than the value of 0.017 in an equilibrium TBL. 
These high values of the mixing length and eddy viscosity 
in the outer parts of the flow far downstream indicate the 
slow response of the Reynolds stresses to changes in the 
mean rate of strain. However, a recovery trend can be 
clearly observed. 

The dissipation length scale put forward by Townsend 
(1961) and cast in somewhat modified form by Bradshaw, 
Ferriss, and Atwell (1967) on the argument that ~uv/k is 
usually well-behaved, is defined as L £ = (-uv) 3/2 /e. It 
can be easily shown that L £ = / , for a flow in energy equi- 
librium. Dissipation obtained by difference of the other 
terms in the balance equation of the turbulent kinetic 
energy was used to calculate L £ . The scatter in the dissi- 
pation length scale, as seen in figure 22, can be attributed 
to the uncertainty of the rate of dissipation. Near the wall 
y < 0.1 6 , L £ increases as kv , as in the case of an equilib- 
rium boundary layer. This observation contrasts with the 
conclusions about Prandtl’s mixing length, but the expla- 
nation is simply that the inner layer is far from energy 
equilibrium. The data show that L £ increases at a higher 
rate for y > 0. 1 5 . In the outer layer, L £ attains peak val- 
ues which are four to five times larger than in an equilib- 
rium TBL. Cutler and Johnston (1989) observed a three- 
fold increase of the same quantity in the recovery region 
behind a fence. The same feature is observed in recovering 


boundary layers perturbed by a concave, destabilizing, 
curvature (Hoffmann, Muck, and Bradshaw (1985); and 
Smits, Eaton, and Bradshaw (1979)) and in the case of 
streamline divergence of developing boundary layers 
(Saddoughi and Joubert (1991)). The spectacular increase 
of the dissipation length scale in the mid-layer is associ- 
ated with the reduced dissipation in the central portion of 
the recovering boundary layer. 


The data indicate that Prandtl mixing length and the dissi- 
pation length parameter gradually recover as the flow 
evolves downstream, but are still far from the equilibrium 
plane boundary layer distributions at the last measuring 
station. Hunt, Spalart, and Mansour (1987) proposed a 
model for the dissipation length 


L e = 


r -u V y 2 (\ dU/dy Vl 

~2 v + As “51/2 

V v ) V y v / 


invoking the local turbulence parameters, in an attempt to 
predict a wider range of wall-bounded complex flows. The 
model underestimates the peak value in mid-layer (not 
shown here). In addition, the model fails to capture the 
decay rate of the dissipation length scale in the stream wise 
direction. 


It appears that P/e , the ratio of the turbulence energy pro- 
duction rate to its dissipation rate, is an important parame- 
ter used for improvement of baseline turbulence models 
for predictions of complex flows (Coakley and Huang 
(1992)). Three different forms of P were determined from 
the data. 



and appropriate ratios are shown in figure 23(a)-(c). The 
structural parameter a = 0.3. It is obvious that P 3 is not an 
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appropriate representation of turbulent kinetic energy pro- 
duction near the wall because the ratio P / £ overshoots 
unity near the wall. It appears that all three forms of P are 
identical in the outer flow region since the production due 
to the normal stresses is negligibly small, while the shear 
stress and kinetic turbulent energy are proportional. The 
ratio P/E dips below unity near the wall, while it over- 
shoots unity in the outer flow region. For large streamwise 
distances, it approaches unity across the entire boundary 
layer as it is expected since the flow approaches equilib- 
rium conditions. 

Conclusions 

The results presented and discussed in the previous sec- 
tions led to the following conclusions about the response 
and recovery of the turbulence structure of the separated 
flow downstream of a backward-facing step. 

The results suggest a two-layer structure of the separated 
boundary layer downstream of the step. The internal layer, 
imbedded in the original boundary layer, grows immedi- 
ately from the step edge while a second layer, the outer 
layer, remains largely unaffected by the disturbance at the 
step, thus resembling the original boundary layer. The 
internal layer behavior is similar to that of a plane mixing 
layer; however, it fails to resemble a mixing layer exactly, 
because of the additional effects of adverse pressure gradi- 
ent, the high turbulence levels on the low-speed side, the 
sheared turbulent “boundary layer” on the high-speed side, 
the curvature, and the short development length before the 
flow reattaches. The boundary between the internal layer 
and the outer region, fed by the original boundary layer, is 
marked by a clearly identifiable kink in the profiles of all 
the measured turbulence quantities. The large structures of 
the two layers appear to have independent dynamics inter- 
acting through the small scales which reside in the inter- 
face between the two layers. Only the presence of the 
small scales between the two layers can support a sudden 


change of gradients detectable in the relevant profiles. Two 
regions are identified in the separated shear layer: the near 
field, x < 2 h, where the mixing-layer-like internal layer 
adjusts to the new boundary conditions, and the far field, 
2h, where vigorous interaction between the internal 
and external layers takes place. The results indicate that 
the two layers merge before the flow reattaches. 

In the recovery region, the new boundary conditions 
imposed on the flow in the reattachment region, u = 0 and 
v = 0, restrict further growth of the large structures of the 
mixing-layer-like flow, and promote a new internal 
boundary layer downstream of reattachment. Based on the 
rate of recovery of the internal boundary layer, two regions 
are identified downstream of reattachment: the near field, 
which ends at x « 20h , where the structure of the internal 
layer attains a quasi-equilibrium state, and the far field for 
x > 20 h approximately. 

The results indicate that the structure near the wall becomes 
wall -dominated when the maximum u rm / u x in the 
external layer decreases to a value of about 3.5. The 
x position where this occurs is a function of two 
parameters, Re h and a perturbation strength, 8/ h . For this 
experiment (Re h = 37,000 and 8/h = 0.82 ) this position 
is x ~ 20 h . If the external perturbation (in the mid layer) 
exceeds this threshold level (which can be also cast in the 
form u /U >3.5 IC,/2) the near- wall structure will be 
strongly altered by the external outer dynamics. It appears 
that, once recovered, the internal layer controls the skin 
friction coefficient. 

Longitudinal turbulent diffusion and production of the 
turbulent energy by the normal stresses (both negligible 
according to the boundary-layer approximation) play a 
significant role in the balance of the turbulent kinetic 
energy in the near field region. Turbulent diffusion and 
dissipation are balanced in the internal boundary layer up 
to about x = 20/i, suggesting that energy equilibrium 
concepts of turbulence modeling are not applicable in this 
region. The importance of the turbulent diffusion in the 
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balance of turbulent kinetic energy and shear stress 
transport equation suggests that an accurate prediction of 
the separated/reattached flow will require a good model for 
the turbulent diffusion terms. In the far field, x > 20 h, the 
internal layer spreads away from the wall, remaining 
self-similar (when scaled on the wall variable) but 
interacting vigorously with the external layer. The 
turbulence quantities in the outer layer decay at a slow rate 
and the external layer still has a memory of the upstream 
disturbance even at the last measuring station of 5\h. 

In the near field (X r <x< 20 h), the nonequilibrium state of 
the flow structure of the internal boundary layer yields a 
linear variation of the mixing length and eddy viscosity 
with y near the wall. The slopes are, however, larger than 
the von Karman constant, K = 0.41, for the mixing 
length and k m x for the eddy viscosity. After the internal 
layer has attained equilibrium (in the far field, x > 20/t), 
the mixing length and eddy viscosity distributions resem- 
ble those of a plane TBL within the developing internal 
layer. Mixing length and eddy viscosity in the external 
layer of the flow exceed the values of a regular TBL by up 
to four times even at the last measuring station. 
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Figure 5. Continued, (c) v 2 . All symbols are same as in figure 5(a). 
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(a) 

Figure 7. Profiles of the internal mixing layer in similarity coordinates, (a) Shear correlation 
coefficient. R uv . All symbols are same as in figure 5(a). 
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(a) 

Figure 8. Profiles of the internal mixing layer in similarity coordinates, (a) Eddy viscosity. All 
symbols are same as in figure 5(a). 
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(a) 


Figure 9. Profiles of the internal mixing layer in similarity coordinates: triple-product transport 
velocities, (a) For turbulent kinetic energy v* =0.6 ^m 
same as in figure 5(a). 
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Figure 10. Continued, (c) Turbulence intensity in the transverse direction. All symbols are same as 
in figure 10(b). 
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Figure 1 1. Development of maximum values in the reattached flow, (a) u^s = / local maximum 

in outer layer, and value aty+ = 20. n~w equil, near wall equilibrium; o~l equif, outer layer 
equilibrium . 

















Figure 12. Concluded, (d) u 
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Figure 16. Development of the internal boundary layer; quadrant distributions of < uv >,• / uv . 
(a) x = -1.05h. 
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Figure 16. Continued, (b) x = 9.87h. All symbols are same as in figure 15(a). 









0.50 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

y/<5 

(b) 

Figure 1 7. Continued (b) x = 9.87h. AH symbols are same as in figure 1 7(a). 


74 








0.01 0.10 1.00 

y/<S 


(a) 

Figure 18. Development of the internal boundary layer budgets of turbulent kinetic energy , 
normalized by u\ /h, with dissipation by difference. Lines denote local equilibrium, 
production = dissipation = u\ / k}\ (a) x = 9.87h. 
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Figure 19. Development of the internal boundary layer: budgets of turbulent shear stress , 
normalized by / 5 , with pressure-strain “redistribution” by difference, (a) x = 9.87h. 
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Figure 20. Development of the internal boundary layer profiles of mixing length and eddy viscosity. 
Dash line represents upstream TBL. (a) Mixing length (sold line presents 0.41 y). 
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(a) 

Figure 23. Development of the internal boundary layer: ratio of the turbulent kinetic energy 
production rate to dissipation rate, (a) P^e. All symbols are same as in figure 20(a). Dash line 
represents equilibrium, production = rate of dissipation. 
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